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o 
o 

O^l ' Exact solvability of two typical examples of the discrete quantum mechanics, i.e. the 

, dynamics of the Meixner-PoUaczek and the continuous Hahn polynomials with full 

D ' parameters, is newly demonstrated both at the Schrodinger and Heisenberg picture 

, levels. A new quasi exactly solvable difference equation is constructed by crossing 

these two dynamics, that is, the quadratic potential function of the continuous Hahn 
polynomial is multiplied by the constant phase factor of the Meixner-PoUaczek type. 
Its ordinary quantum mechanical counterpart, if exists, does not seem to be known. 



1 Introduction 



I As shown recently, Quasi Exact Solvability (QES) is very closely related to exact solvabil- 

ity [22\ \T7\ |2T] . If all the eigenvalues of a quantum mechanical system are known together 
with the corresponding eigenfunctions, the system is exactly solvable in the Schrodinger 
I picture. In contrast, a system is QES if only a finite number (usually the lowest lying 

ones) of exact eigenvalues and eigenfunctions are known [26^ [25} I24j . Among various char- 
■ acterisation of quasi exact solvability \26\ [25| [2^ [T2l [8], the existence of the invariant 

^ i polynomial subspace is conceptually simple. The method to obtain a QES system, advo- 

CN I cated by the present author [22l [TTJ [21] , by deforming an exactly solvable system with an 

addition/multiplication of a higher order interaction term together with a compensation 
(3 ■ term, exemplifies the structure of the invariant polynomial subspace rather clearly through 

the action of the similarity transformed Hamiltonian H (j3.5p - (j3.7p in terms of the pseudo 
groundstate wavefunction This method was applied to the exactly solvable ordinary 
quantum mechanics [21] of one degree of freedom and multi-particle system of Calogero- 
I Sutherland type [4l[23]. Recently new QES difference equations of one degree of freedom 

[22j and multi-particle systems [17J are obtained by the application of the same method to 
the discrete quantum mechanics [14] for the Askey-scheme of hypergeometric orthogonal 
polynomials [B [H] and for the Ruijsenaars-Schneider-van Diejen systems \20\ I27j. 

Two of the exactly solvable discrete quantum mechanics discussed in [HI [HI [22] , the 
Meixner-Pollaczek and the continuous Hahn polynomials, are of special types in the sense 
that their parameters subset of the allowed ones. 

The purpose of the present paper is three-fold. Firstly, to demonstrate the exact solv- 
ability of the full dynamics of the Meixner-Pollaczek and the continuous Hahn polynomials 
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in the Schrodinger picture through shape invariance [7l[T3j. The exact Heisenberg oper- 
ator solutions are also constructed through the closure relations (j2.12p . (j2.26p . (j2.38p . 
The structure of the invariant polynomial subspace is shown explicitly by the action of 
H for each degree monomial (j2.22p . (|2.33p . Secondly, to obtain a new QES difference 
equation by crossing the above mentioned exactly solvable dynamics. The new system has 
the quadratic potential with two complex parameters ()2.29p coming from the continuous 
Hahn polynomial and a constant multiplicative phase factor e~'^^ (I2.17P coming from the 
Meixner-Pollaczek polynomial. Thirdly, to give comments on exact Heisenberg operator 
solutions. The third part is closely related to the presentation in NEEDS 2007 Workshop 
by the present author, "Heisenberg operator solutions for the Calogero systems" [16j. 

This paper is organised as follows. In section two the exact solvability of the full 
dynamics of the Meixner-Pollaczek and continuous Hahn polynomials is demonstrated 
after brief review of the general setting of the discrete quantum mechanics appropriate for 
the Askey-scheme of hypergeometric orthogonal polynomials. Section three is devoted to 
the new QES difference equation obtained by crossing the dynamics of the full Meixner- 
Pollaczek and continuous Hahn polynomials. Section four is for the comments on the 
exact Heisenberg operator solutions. Their dynamical roles, algebraic interpretation and 
the connection to a ^quantum Liouville theorem^ are explained. 

2 Hamiltonian Formulation for Dynamics of Hypergeomet- 
ric Orthogonal Polynomials 

It is well known that the classical orthogonal polynomials, the Hermite, Laguerre and Ja- 
cobi polynomials with various degenerations (Gegenbauer, Legendre, etc) constitute the 
eigenfunctions of exactly solvable quantum mechanics, for example, the harmonic oscil- 
lator without /with the centrifugal potential, the Poschl- Teller potential etc. Thus it is 
quite natural to expect that the Askey-scheme of hypergeometric orthogonal polynomi- 
als together with their g-analogues, which are generalisation/deformation of the classical 
orthogonal polynomials, also constitute the eigenfunctions of certain quantum mechanics- 
like systems, so that the orthogonality has proper explanation/interpretation. In 'discrete' 
quantum mechanics [TJj, a Hamiltonian formulation was introduced for the dynamics of 
several typical examples of the Askey-scheme of hypergeometric orthogonal polynomials. 
Since these polynomials obey difference equations instead of differential equations, the 
Hamiltonians contain the momentum operators in the exponentiated forms in contrast to 
the second order polynomials in ordinary quantum mechanics. These examples of dis- 
crete quantum mechanics are exactly solvable in the Schrodinger picture due to the shape 
invariance properties [3 [H] and their exact Heisenberg solutions are given in [15j . 

In this section we discuss two examples, the Meixner-Pollaczek polynomial and the 
continuous Hahn polynomial, in their full generality. In our previous work on discrete 
quantum mechanics \14:\ I15j . only the special case of the Meixner-Pollaczek polynomial 
with the phase angle ^ = 7r/2 and the special case of the continuous Hahn polynomial 
with two real parameters ai and 02 are discussed, partly because these special cases of the 
two polynomials appear in several other dynamical contexts [Sl [U [2] and, in particular, 
they appear as describing the equilibrium positions [19l[13l[l3j of the classical Ruijsenaars- 
Schneider van Diejen systems pOt [27]. 
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We sliow tliat these two polynomials in their full generality, that is, with a general 
phase angle (p for the Meixner-Pollaczek polynomial and with two complex parameters oi 
and a2 for the continuous Hahn polynomial, are exactly solvable in the Schrodinger as 
well as in the Heisenberg picture. Later in section [3] we show that a new quasi exactly 
solvable system is obtained by crossing these general Meixner-Pollaczek and continuous 
Hahn polynomials, that is, by multiplying the extra phase factor to the potential of the 
general continuous Hahn polynomial. The resulting system is no longer exactly solvable 
but it is quasi exactly solvable by adding a compensation term. 



2.1 General Setting 

Here we recapitulate the basic notation and rudimentary facts of discrete quantum me- 
chanics of one degree of freedom. For details we refer to |14l I18j . The Hamiltonian of 
a discrete quantum mechanical system of one degree of freedom to be discussed in this 
paper has the following general structure 

n = VV{x) e^^V{xf + ^/V{xY e"P^V{x) - V{x) - Vix)* (2.1) 
= \/V(x) e"*^- ^V(x)* + ^V{xY e+*^^ -JVix) - V(x) - V{xf, (2.2) 

in which the potential function V(x) = V{x ; A) depends, in general, on a set of parameters 
A. The exponentiated momentum operators cause a finite shift of the wavefunction in the 
imaginary direction: e^'^''(/>(x) = (/)(xibi). As in the supersymmetric quantum mechanics 
[9l [5] , the Hamiltonian is always factorised 

n = A^A, (2.3) 
A^ = ^/V^e-^^^-^/V{xrei^\ A = e-i^^^/V{xr-ei^^^/V{x), (2.4) 

which shows the (formal) hermiticity and positive semi-definiteness of the Hamiltonian. 
See the discussion in §4 of [22j for detailed realisation of hermiticity. The groundstate 
wavefunction (j)Q{x) is annihilated by the A operator 

A(l)o{x) = =^ TiMx) = 0, (2.5) 

which can be chosen real (poix) G The eigenfunctions of the Hamiltonian (pnix) = 
(f)n{x ; A) have the following general structure: 

nMx)=£nMx) (n = 0,1,2,...), = ^0 <^i <<?2 < ••• , (2.6) 
</.„(x;A) = (/.o(x;A)P„(r/(x);A). (2.7) 

Here Pn is a polynomial in r]{x), which is called a sinusoidal coordinate [15]. The orthogo- 
nality theorem for the eigenfunctions belonging to different eigenvalues implies that {Pn} 
are orthogonal polynomials with respect to the weightfunction 0o(^)- 



/ 



^o(x ; X)Pnirj{x) ; X)*Pm{rjix) ^)dx oc 6nm- (2.^ 
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Shape Invariance If the reversed order Hamiltonian AA^ has the same form as the 
original Hamihonian A'^A 

^(A)^(A)t = ^(A + d)^A{X + S)+ £i{X), (2.9) 

the system is called shape invariant [Tj [14]. Here S denotes the shift of the parameters 
and an additive constant <?i(A) is to be identified as the energy of the first excited level. 
Combined with the basic fact of the super symmetric quantum mechanics that the two 
Hamiltonians A'^A and AA^ are iso-spectral (except for the groundstate) , shape invari- 
ance determines the entire energy spectrum and the excited state eigenfunctions from the 
groundstate wavefunction: 

£nW = Y.£i{X + sS), (2.10) 

s=0 

4>n{x ; A) oc A(A)U(A + S)^A{X + 26)^ ■■■A{X + {n- l)d)Uo{x ; A + nd). (2.11) 
This establishes the exact solvability in the Schrodinger picture. 



Heisenberg Operator Solution The sinusoidal coordinate 'r]{x) has a remarkable prop- 
erty [15j that the multiple commutators with the Hamiltonian can be reduced to r]{x) itself 
and the first commutator [TC, 7]] through the closure relation 

[n,[n,7]]] = 7]Ro{n) + [n,v]Ri{n) + R-i{n). (2.12) 

Here Rq{H) and R-i{'H) are in general quadratic polynomials in TL, whereas RiiTL) is lin- 
ear in 7i. This leads to the exact Heisenberg operator solution for the sinusoidal coordinate 
r/(x): 

e^*^??(x)e-**^ = aWe^°+(^)* + a^-^e^"-^^)* - i?_i(H)/i?o(W), (2.13) 

a±{n) = \ {Ri{n) ± ^jRiiny + mo{n) ) , (2.14) 
„(±) (±[n,v{x)] T {7]{x) + R^i{n)/Ro{n))a^{n)^ / {a+{n)-a^{n)). (2.15) 

The entire spectrum {£n} can also be determined from (j2.13p by starting from £o = jT5] . 
as done by Heisenberg and Pauli for the harmonic oscillator and the hydrogen atom. The 
positive and negative energy parts a^^^ of the Heisenberg operator solution e**^r;(a;)e~**^ 
are the annihilation-creation operators: 

Q(+)t = a(+)(/)„(x) oc (j)n+i{x), a("Vn(a;) oc (j)n~i{x). (2-16) 



The general theory of exact Heisenberg operator solutions for exactly solvable multi- 
particle systems is yet to be constructed. For the special case of the Calogero systems 
[U [To] , the totality of Heisenberg operators are derived for any root systems [K] . For the 
classical root systems A, BC and -D, the number of particles can be as large as wanted. 
See section [4] for comments on exact Heisenberg operator solutions in general. 
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2.2 Meixner-Pollaczek polynomial 

Tlie potential function V{x) for the Meixner-Pollaczek polynomial is linear in x: 

V{x;\)=e-'f^{a + ix), X = a, (2.17) 
0<aGM, 0gR, (3=^(j)-^, 0<(/><7r. (2.18) 

The special case discussed in [HI |6l[2l[3] is/3 = 0or(/) = 7r/2. The groundstate wavefunc- 
tion (pQ, as annihilated by the operator A, A(pQ = 0, is given by 

; a) = e^^|r(a + ix)\ = e'^^-^^'' y^T{a + ix)T{a - ix). (2.19) 
The similarity transformed Hamiltonian H in terms of the groundstate wavefunction ^o, 

n = (j)-^ono^Q = v{x) (e-'^^ - 1) + v{xy (e*^- - 1) (2.20) 



(a + ix) e-'^ (e-'^^ - l) + {a - ix) e'^ (e'^^ - 1 ) (2.21 



acts on the polynomial part of the wavefunction. It is obvious that 7i maps a polynomial 
into another and it is easy to verify 

Hx^ = 2ncos P x"' + lower order terms, n € (2.22) 

Thus we can find a degree n polynomial eigenfunction Pn{x) of the similarity transformed 
Hamiltonian H. 

nPn{x) = £nPn{x), £n = 2n COS P = 2n siu (j), n = 0,1,2,..., (2.23) 

which is called the Meixner-Pollaczek polynomial [TT]. It is expressed in terms of the 
hypergeometric series 

pW(^;^) = Mie-*2i7,(-^'« + ^^ l-e-2.<^), (2.24) 



n! V 2a 

in which (b)n is the standard Pochhammer's symbol 



ib)n = ll{b + k-l)=b{b + l)---{b + n- 1). 



def 

k=l 



Shape invariance is also easy to verify: 

A{x;a)A{x;a)^ = A{x;a + l)^A{x;a + l)+£i{X), fi(A) = 2sin</). (2.25) 

The parameter a is increased by ^ whereas the new parameter (f> is invariant. Since 
£1 is independent of the shifted parameter a, it is trivial to obtain the linear spectrum 
£n = 2nsm(j), which is the same as (12.23p . 

The sinusoidal coordinate is ry(a;) = x. The closure relation (I2.12P reads simply 

[n,[n,x]] = X A sin'^ + 2 cos (t>n + 2a sin 2(1), a±{n) = ±2 sin 4>, (2.26) 
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indicating that x undergoes a sinusoidal motion with the frequency 2 sin (p. The Heisenberg 
operator solution is 



jtH^-it'H_^ rn.„;^^i , .rov ^iSin[2tsin( 



^ttn ^ ^-ttn ^ ^ pQgj2t sin (p] + i[n, x 



2 sin 

cos 



2 sin^ , 

The annihilation-creation operators are: 



{n + 2asin(/>)(cos[2tsin(/>] - 1). (2.27) 



a^^) =±[n,x]/{4sm^) + -\x+ ^°^f (H + 2a sin (/.) I . (2.28) 

2 [ 2 sin (f) J 

Obviously the expressions (I2.26p - ()2.28p are drastically simplified for the special case of 
= 7r/2 which were discussed in previous work \14:\ ITS] . 

2.3 Continuous Hahn polynomial 

The potential function V{x) for the continuous Hahn polynomial is quadratic in x: 

V{x ] X) {ai+ix){a2+ix), A = (ai, 02), oi, 02 G C, Re(ai) > 0, Re(a2) > 0. (2.29) 

The special case discussed in [HI El El |3] is for real oi and a2- The groundstate wavefunc- 
tion as annihilated by the operator A, AcpQ = 0, is given by 



; A) =^ Y'r(ai + ix)T{a2 + ix)T{al - ix)r(a^ - ix). (2.30) 
The similarity transformed Hamiltonian H in terms of the groundstate wavefunction (pQ, 



n = <p-^ono(t>^ = V{x) (e-'^^ - 1) + V{xY (e'^- - 1) 



(2.31) 



(ai + ix){a2 + ix) ( e~'^- - 1 ) + K - ix){al - ix) ( e*^- - 1 ) (2.32) 



acts on the polynomial part of the wavefunction. It is obvious that 7i maps a polynomial 
into another and it is easy to verify 

= n(n + ai + 0*1 + a2 + — l)x"' + lower order terms, n € Z+. (2.33) 

Thus we can find a degree n polynomial eigenfunction Pn{x) of the similarity transformed 
Hamiltonian 7i 

nPn{x) = £nPn{x), f„ = n(n + oi + + 02 + - 1), n = 0,1,2,..., (2.34) 

which is called the continuous Hahn polynomial [TT]. It is expressed in terms of the 
hypergeometric series 



p„(x; 01,02, 01,02) (2.35) 

iV 



.„ (oi + a^)„(oi + a2)n ( -n, n + oi + 02 + oj; + 02 - 1, oi + 
n! V oi + ai, oi + Og 
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Shape invariance is also easy to verify: 

A{x;ai,a2)A{x;ai,a2)^ = A{x;ai + ^,02 + ^)'^A{x;ai + ^,02 + ^) +<?i (01,02), 

(2.36) 

(5=^(i,i), <Si(ai, 02) = fei, 61 =^ ai + 02 + a* + 02 = 2Re(oi + 02). (2.37) 

Here we have introduced an abbreviation bi for convenience. The parameter oi and 02 are 
increased by ^. Since is hnearly dependent on the shifted parameters ai, and 02, it is 
trivial to obtain the quadratic spectrum = ^{71 + oi + 02 + + 02 — 1) = n(n + 61 — 1), 
which is the same as (|2.34p . 

The sinusoidal coordinate is r]{x) = x. The closure relation (I2.12P reads simply 

[n, [n, x]] =xim + hih - 2)) + 2[n, x] + 62W + 63(^1 - 2), (2.38) 

in which abbreviations 62 '= 2Im(ai + 02) and 63 '= 2Im(aia2) are used. The frequencies 
a±{Ti.) are 

a±{n) = 1 ± 2Vn', n' = n + {h - n'<Pn = {n + {h - i)/2fcpn- (2.39) 

The Heisenberg operator solution reads 

^UH^^-an -a-(7^)e-+(^)^ + a4^)e--(^)* ^ e-+(>^)* - e--(^)* 



+ 



+ 63(61 - 2) /-a_(W)e^"+(^)* + a+(W)e^"-(^)* 
4(W + 6i(6i-2)) I iTw^ 



(2.40) 



The annihilation and creation operators are: 

„/(±)def^(±)4^ 



4(?^ + 6l(6l-2)), 

Obviously the expressions (j2.38p - (|2.4ip become drastically simple for the special case of 
62 = 63 = 0, which were discussed in previous works |141 115j . 



3 New QES Difference Equation 

Here we will discuss the discrete quantum mechanics obtained by crossing the Meixner- 
Pollaczek and the continuous Hahn polynomials, that is, with the quadratic potential 
function of the continuous Hahn polynomial ()2.29p multiplied by a constant phase factor 
g-«/3 Qf i^j^g Meixner-Pollaczek type. As vaguely expected, the exact solvability is not 
realised. We will show, instead, that the system is quasi exactly solvable by adding a 
compensation term which is linear in x: 

n .JV{x)e-'^^^/V^* ^ ^V(^ e^'^^^V{^-V(x) - V(xy + a^x (3.1) 

= ylU + a^x, =^ -2A^ sin /?, 7W G Z+, (3.2) 

y(x) =^ (ai +ia;)(o2 + ix)e~*^, ai,a2GC, Re(ai) > 0, Re(a2) > 0. (3.3) 
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It should be noted that the Hamiltonian is no longer positive semi-definite but the her- 
miticity is preserved. The main part, that is without the compensation term, is factorised 
as before (12. 4p . The zero mode of the A operator 

A(po = 0^ Mx) = e^'^yr(ai + ix)T{a2 + ix)T{al - ix)T{a*2 - ix), (3.4) 

is no longer the groundstate wavefunction. It is called the pseudo groundstate wavefunction 
[22]. 

The similarity transformed Hamiltonian 7i in terms of the pseudo groundstate wave- 
function (j)Q, 

n = (j)^^ on 0^0 = V{x) (e-*^- - l) + V{xy (e*^- - l) + aMX (3.5) 

= (ai + ix){a2 + ix)e-'l^ (e"*^- - l) + {a[ - ix){al - ix)e''^ (e^^^ - l) 

-2Msmpx, (3.6) 

acts on the polynomial part of the wavefunction. It is obvious that 7i maps a polynomial 
into another and it is easy to verify 

= 2{-M + n) sin P + lower order terms, n G Z_^. (3.7) 

This means that the system is not exactly solvable without the compensation term, but 
it is quasi exactly solvable, since 7i has an invariant polynomial subspace of degree Ai: 

HVmQVm, (3.8) 

=^Span[l,x,x2,...,2;^] , dimVA^ = 7W + 1. (3.9) 

The Hamiltonian 7i ()3.2p is obviously hermitian (self-adjoint) and all the eigenvalues are 
real and eigenfunctions can be chosen real. We can obtain a finite number (A^ + 1) of exact 
eigenvalues and eigenfunctions for each given A4 . The oscillation theorem linking the num- 
ber of eigenvalues (from the groundstate) to the zeros of eigenfunctions does not hold in the 
difference equations. The square integrability of all the eigenfunctions (jp'{x)dx < oo 
is obvious. See [22j for other examples of quasi exactly solvable difference equations of one 
degree of freedom and [T7] of many degrees of freedom. 

It is easy to demonstrate that multiplying an extra constant phase factor e~*^ to the 
other exactly solvable potential functions 



def {a,+^x){a2 + ^x){a, + ^x) continuous dual Hahn, (3.10) 

def(ai+.x)(a, + .x)(a3 + ix)(a4 + ix)^ Wilson, (3.11) 
2ix{2ix + 1) 

V{x) g - ''''^^] - "y,^/^ - - "^"^ , z = e«, Askey-Wilson, (3.12) 
(1 — z'^){\ — qz-^) 



does not provide either exactly solvable or quasi exactly solvable dynamical systems. The 
situation is the same for various restrictions of the Askey-Wilson polynomial. 
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It is also easy to see tliat the quasi exact solvability of the systems discussed in [22] 
with the potentials 

V{x) =^ (ai + ix){a2 + ix){a3 + ix), (3.13) 

V{x) =^ (ai + ix){a2 + ix){a3, + ix){a4 + ix), (3-14) 

. dei {ai+ix){a2+ix){a3+ix){a4 + ix){a5+ix) /o r\ 

^("^ = 2^x(2^x + 1) ' ^''''^ 

X def (qi + ix){a2 + ix){a3 + ix){ai + ix){a^ + ix)(a6 + ia:) ^'^^p^ 

~ 2ix{2ix + 1) ' ^ ■ 
, . def (1 - aiz)(l - a2z){l - azz){l - aiz){l - a^z) _ 
(1 — 2^)(1 — qz'^) 

is destroyed if a constant phase factor e~*^ is multiplied. 



4 Comments on Exact Heisenberg Operator Solutions 

Let us start with a rather naive question; "What can we learn more from the exact Heisen- 
berg operator solutions when we know the complete spectrum and the corresponding eigen- 
functions?" A small digression on the well known relationship between the Schrodinger 
and Heisenberg pictures would be useful. Suppose we have a complete set of solutions of 
the Shrodinger equation 

For any observable A, one can construct a, (usually infinite) matrix yl, Ayitti — (^nl^l^Tn); 
satisfying the Heisenberg equation of motion 

Obviously such an exact Heisenberg operator solution does not teach us anything more. 

But for a special choice of the observables, called the 'sinusoidal coordinates' {r]j}, 
j = 1, . . . ,r, with r being the degree of freedom, the operators 

{e*%,e-^^*}, i = l,...,r, 

can be expressed explicitly in terms of the fundamental operators {r]j}, Ti. and a finite 
number of multiple commutators of {rij} with the Hamiltonian [7i, [7i, [■ ■ ■ ,?7j]..]. These 
are the Heisenberg operator solutions found by Odake-Sasaki for a wide class of exactly 
solvable degree one quantum mechanics including the discrete ones [15] and for typical 
multi-particle dynamics of Calogero type for any root system [16]. It should be stressed 
that the existence of sinusoidal coordinates is not guaranteed at all. There are several 
exactly solvable degree one quantum mechanical systems for which our construction of 
the Heisenberg operator solutions does not apply. Various reduced Kepler problems and 
the Rosen-Morse potentials are the typical examples. See |15] for more details. For multi- 
particle systems, the exact Heisenberg operator solutions are known [16] only for the 
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Calogero systems for any root system [HHQ]. There are other well-known exactly solvable 
multi-particle systems; the Sutherland systems [23] and the Ruijsenaars-Schneider-van 
Diejen systems [201 127] , The name 'sinusoidal' implies that they all undergo sinusoidal 
motion but not harmonic. In classical mechanics terms, the frequencies depend on the 
initial conditions. 

From the analysis point of view, the sinusoidal coordinates generate the polynomial 
eigenfunction {Pn}, (j^n = <t'oPn (0o is the groundstate wavefunction) . In other words 
{Pn} are orthogonal polynomials in {rjj}. The exact Heisenberg operator solutions for 
{r/j} provide the complete set of multi- variable generalisation of the three term recurrence 
relations, which characterise orthogonal polynomials in one variable. 

As stressed in 116] . the positive and negative frequency parts of the Heisenberg 
operator solutions are the sets of annihilation- creation operators. They generate the entire 
eigenfunctions algebraically, and thus form a dynamical symmetry algebra together with 
the Hamiltonian and possibly with the higher conserved quantities (Hamiltonians). The 
structure of these dynamical symmetry algebras is identified only for a few special cases 
of degree one, for example, su(l,l). It is a good challenge to identify the dynamical 
symmetry algebra and its irreducible representations for each known exact Heisenberg 
operator solution, single and multi-degree of freedom. From the algebra point of view, 
the three term recurrence relations for single variable orthogonal polynomials correspond 
to the Clebsch-Gordan decomposition rules for rank one algebras. The multi-particle 
version would simply correspond to the higher rank counterparts of the Clebsch-Gordan 
decomposition rules. 

From a more basic dynamics point of view, one could consider the exact Heisenberg 
operator solutions and the associated annihilation-creation operators as an explicit but 
partial realisation of ^quantum Liouville theorem^ The classical Liouville theorem asserts 
that one can construct by quadrature only from the complete set of involutive conserved 
quantities the generating function of a canonical transformation which brings the system 
to the action-angle form. In contrast, the usual formulation of quantum Liouville theorem 
does not say anything about the second half; quantum mechanical counterpart of 'bringing 
to the action-angle form'. The complete set of the creation-annihilation operators play 
the corresponding role; 'generating the entire eigenfunctions from the groundstate wave- 
function'. If such generated eigenstates were the simultaneous eigenstates of the complete 
set of involutive conserved quantities, one could say that the quantum Liouville theorem 
is fully realised. It seems that there is still some way to go for that goal. 
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